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The Chapman-Enskog method is generalized for accounting the effect of kinetic modes on hydrodynamic evolu- 
tion. Hydrodynamic states of phonon system of insulators have been studied in a small drift velocity approxima- 
tion. For simplicity, the investigation was carried out for crystals of the cubic class symmetry. It has been found 
that in phonon hydrodynamics, local equilibrium is violated even in the approximation linear in velocity. This 
is due to the absence of phonon momentum conservation law that leads to a drift velocity relaxation. Phonon 
hydrodynamic equations which take dissipative processes into account have been obtained. The results were 
compared with the standard theory based on the local equilibrium validity. Integral equations have been ob- 
tained for calculating the objects of the theory (including viscosity and heat conductivity). It has been shown that 
in low temperature limit, these equations are solvable by iterations. Steady states of the system have been con- 
sidered and an expression for steady state heat conductivity has been obtained. It coincides with the famous 
result by Akhiezer in the leading low temperature approximation. It has been established that temperature 
distribution in the steady state of insulator satisfies a condition of heat source absence. 
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1. Introduction 

A lot of papers are devoted to investigations of the phonon system of insulators. We mean hydrody- 
namic states described by densities of additive conserved values or their functions (for a system with 
broken symmetry, order parameters are added). The description of hydrodynamic processes in insula- 
tors is still more complicated due to the presence of Umklapp processes in which the momentum is not 
conserved. A pioneering research into the role of such processes in the theory of thermal conductivity of 
insulators belongs to Peierls |1]. Significant contribution to the study of the problem has been made by 
Akhiezer [2j who calculated the thermal conductivity of dielectric crystals in steady states at low tempe- 
ratures. 

In phonon hydrodynamics, local temperature T{x, t) and drift velocity u n {x, t) introduced with the 
help of the local equilibrium distribution, are used instead of the densities of energy e(x, t) and mo- 
mentum n n {x, t) (see, for example, books and reviews (sj-Ql containing the standard theory). The 
mentioned works substantiated that the velocity decays with time (slowly at low temperatures) but this 
relaxation process is not considered to be a manifestation of the kinetic mode of a system. During the 
recent years, great attention has been paid to the investigation of the effect of kinetic modes on the sys- 
tem dynamics. Especially the role of the Lviv school on statistical physics (Mryglod, Bryk, Tokarchuk, 
Omelyan, et al.) should be mentioned. See in this regard, for example, the paper (8[] devoted to the effect 
of relaxation on hydrodynamic processes in a two-component system, and the review @] on the general 
theory. Similar studies of dielectrics are unknown for the author and it was the motive of the investiga- 
tion. In the investigation we observed that the local equilibrium distribution at u n {x, t) + does not give 
the leading approximation for the nonequilibrium distribution function. Similar result concerning the 
violation of the local equilibrium was obtained in our paper fioll on the polaron kinetics due to the pres- 
ence of velocity and temperature relaxation in the electron-phonon system. Such a result has not been 
obtained in (J since the dynamics of the system was studied in terms of component energies. 
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Researchers usually put the Chapman-Enskog method in the basis of hydrodynamic state considera- 
tion. In the general formulation of this method, the particle distribution function f ap (x, t) (a is the num- 
ber of species of particles, p n is momentum) is considered to be a functional f ap (x,£ (r)) of parameters 
£ p (x, t) describing the hydrodynamic state (see, for example, I'll]). Such state is observed after the time 
To that depends on the initial state of the system and the problem to be solved 

f ap (x,t) -W*,f(f)). (1) 

f»T 

Relation (1) is referred to as the functional hypothesis and is a natural assumption allowing one to obtain 
a closed set of equations for parameters £ p (x, t) and to express all the observed values through these 
parameters. In the known paper |2], Akhiezer actually made the same assumption concerning the struc- 
ture of the nonequilibrium phonon distribution function in his consideration of stationary states of the 
phonon system of insulators. 

A general study of hydrodynamic states of the phonon system of dielectric crystal is conducted in the 
standard theory \^[f^ using the Chapman-Enskog method. Herewith, the coincidence of the distribution 
function f ap (x,£(r)) in the zero approximation in the gradients of hydrodynamic variables f^„(x,£(?)) 
and the Planck distribution with the drift velocity u n {x, t) and temperature T{x, t) 

n ap (£{x,t)), n ap & = [e [£ -p-P" u - )IT (2) 

is substantiated [^(x, t): £o(x, t) = T{x, t), £„{x, t) = u n {x, £)]. The relation f<f p (x,if) = n a p[^[x)) means 
that the local equilibrium is valid for the system. At low temperatures, this result is grounded with great 
accuracy by the fact that the number of the Umklapp collisions is small compared to the number of the 
normal ones. In the present paper, the zero approximation in the gradients f^pix,^) for the distribution 
function f ap (x,f) is calculated not only for low temperatures. This is done with the appropriate gener- 
alization of the Chapman-Enskog method and, in the case of small drift velocity u n , it is shown that 
C P (x,0 * n ap &x)) but C P (x,0\u=o = n ap {${x)) (n ap (0 = [e e "P IT - 1]" 1 is the local Planck distribution 
without velocity). In all the mentioned works , the drift velocity of phonons u n is also considered to 
be a small value but our consideration yields nonlinear hydrodynamics for a small drift velocity [it con- 
tains nonlinear terms in u n {x, t)]. The established violation of the local equilibrium has principal physical 
meaning although at low temperatures it is small. 

In the present paper, insulators of a cubic symmetry classes are considered. Even in this case, some 
crystal quantities have a complicated structure (for example, tensors of the fourth rank). Using a small 
drift velocity approximation, let us avoid a complicated symmetry research. Among other things, terms 
in estimations have a complicated tensor structure in all our series expansions in the velocity. 

In some papers, devoted to phonon hydrodynamics issues, crystalline solids are considered in an 
isotropic model . This leads to some simplification of the equations obtained in the present paper 

that will be analyzed elsewhere. A version of a r-approximation for phonon collision integral proposed 
in fl2ll allows one to simplify the equations of the theory too. Such an approximation is not used in the 
present paper. Moreover, a T-approximation approach can be improved by the method of orthogonal 
polynomials developed in |7|]. 

The outline of the paper is as follows. In section[2] the basic equations of the theory are presented. In 
section[3] the basic equations of the theory are solved in a double perturbation theory in small gradients 
of hydrodynamic variables and small velocity. As a result, expressions for the fluxes of energy and mo- 
mentum of the phonons with allowance for dissipative processes are obtained and equations of phonon 
hydrodynamics are derived. Section|4]discusses the relation of the results with the standard phonon hy- 
drodynamics based on the validity of the local equilibrium. In section[5] an iteration procedure is built for 
solving integral equations of the theory in the approximation of low temperatures. In the final section[6] 
the steady states of the system are considered. 
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2. Basic equations of the theory 

A kinetic equation for phonons is put 

df ap {x, t) _ d£ ap df ap {x, t) 
dt dp n dx n 



+ Iap1Kx.fi) (3) 



in the basis of the consideration. The main contribution to the collision integral corresponds to three- 
phonon processes 

lap® = £ f d 3 p 1 d 3 p2d 3 p 3 mi2,3)\ 2 S(E 1 +£ 2 -E3)Y. S 1Pi + P2-P3-flb n ) 

[abbreviations of the type £ ai pi = f a ,p, = h> < &{tt\P\',a2P2,a , iP , i) = $(12,3) are used; a, p n are polar- 
ization and quasi-momentum of a phonon; the integrals are taken over the basic cell of a reciprocal lattice 
B; for simplicity, the term "momentum" is used instead of "quasi-momentum"]. The collision integral is 
divided into the contribution of normal processes and Umklapp ones 7 ap (f) = /^(f) + ^ p (f) in the usual 
manner [in J4) summand with n = gives I^ p (f}]. The hydrodynamic equations are a consequence of the 
phonon energy conservation law and the law of system momentum change 

de{x,t) dqi{x,t) dn n {x,t) dt„i{x,t) 

= 1 = vR n {x,t). (5) 

dt dxi dt dxi 

Densities of energy and momentum e{x, t), n n [x, t) of the phonon system, corresponding flux densities 
qi{x, t), t n i[x, t), and that of a frictional force R n {x, t) are given by formulas 

£(x, t) = £ l dT p £ ap f ap (x, t), n n {x, t)s^ &x p p n i ap {x, t), 

a J a J 

B B 

dT p £ ap —^f ap (x, t), t n i{x, t) = I fapPn-; —fapix, t), 
"Pi a J "Pi 



R n {x,t) = £ I dT p p n I% p (f(x,t)), 

a J 



B 
43, 



Temperature T[x, t) and phonon drift velocity u n {x, t) [variables ( p {x, t)] are used in this paper in- 
stead of densities e(x, t), n n (x, t) as independent variables in hydrodynamics. In the standard theory 
ma , the densities are expressed through £ p (x, t) with the definitions 

£(«D=Z / dT p E ap n ap {0, n n {^)=Y J I dT p p n n ap {£,) (7) 

a J a J 

B B 

where n ap [^) is the Planck distribution 0. Below, we use the details connected with these relations only 
while comparing the presented theory with the standard theory according to the above remarks of ours 
on the violation of the local equilibrium for u n + 0. 

Investigating fundamental issues of phonon kinetics, we restrict our consideration to the crystals of 
cubic symmetry classes O, O/j, I), for simplicity. In the present paper, the phonon drift velocity u n is 
assumed to be small compared to the second sound velocity (let A be a corresponding small parameter). 
In this case, the energy and momentum densities of the phonon gas can be expanded in powers of u n as 
it follows 

£{0=£{T,u)=£{T) + a{T)u 2 + O{A, 4 ), n n {0 = n n {T, u) = a{T)u n + 0(A 3 ) (8) 
where functions £{T), a{T), a{T) are considered to be known. In the standard theory 

ma, 

they are 

calculated from 0. The terms given in (8) by estimations have a complicated tensor structure even for 
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the mentioned symmetry classes. Their complete investigation can be made using the ideas developed, 
for example, in 1 1311 . 

Hydrodynamic equations for the variables S, p {x, t) have the structure 



dt 



M^[xM(t))), 



d s { p (x,t) 



(9) 



and, as usual, the gradients of these variables are small [in (9) g = II L, where I is a mean free path of 
particles, L is a characteristic size of inhomogeneities in the system]. Functions M p (x,f) can be expressed 
through the distribution function f ap (x, f) with the help of relations ([5), (6), and (|8). According to iQJ and 
(9), the functional f ap (x, <f) satisfies the integro-differential equation 



E f d 3 * 

a J 



and additional conditions 



,Sf ap (x,0 



de ap df ap (x,0 
dp n dx n 



+ I ap {f(x,0) 



(10) 



Y, \ dT p e ap f ap {x,0 = £(f(x)), Y, I dT p p B f op (x,0 =7rn(fM) (11) 



that define the phonon gas temperature and drift velocity together with (8). 



3. Phonon hydrodynamics 

Equations O, (|ll) are solved by us in the double perturbation theory in gradients (the small param- 
eter g) and in the drift velocity (the small parameter A) 



fa P (x, o = C + 4 P + o (g 2 ) , 4 P = 4 P 0) + 4 P X) + o(g s x 2 y, 47 = «* P ■ 



(12) 



The Planck distribution without velocity n ap is the main contribution to the phonon distribution func- 
tion. Our calculation is based on the estimations 



T~X 



o„o 



u„ ~ A 



i „o 



dT 



aV. 



dx„ 

The simplest and the most important contributions to the theory have the following structure 
4V = ^ap[ 1 + ^ap)A an (p)u n , 4°p 2) = »ap(l + «ap)fia»z(p)"7i"Z. 

iap — n apy L ^ ,l apj K -'an\H' q^, > L ap ~ n ap y 1 ^ ,l apj 
These formulas contain functions A an {p), C a n(p), D an i[p) that satisfy equations 

-A an {p) = T / dV^^.p'J^.lp'), <p»Am(p)> = 3<r; 



(13) 



DanlW^— + E an i{p)- — U/ 
OX/ ox„ 



(14) 



(15) 



Aan(p)^fvT + ^y^2£a P ^ L = Y f f P* K a a< [p, p')C a > n{p') , (p n C an {p)) = 0; (16) 

« a' B 

-^2£a p Snl + -D an l{p) - A an {p) = £ f d VW [p, p')D a 'nl [p') , 

cT z r a dpi ^TJ 



(£a p D an „{p)) = 0. 



(17) 
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Here, heat capacity c, pressure p of the phonon gas and a special average (g ap ) 
de ( £ l p ) 1 f , de ap 

c - 



dT T 2 



i c o£ a p o ^ r o i o \ 

p=-ldT p p n - n ap ; (gap) = 2^\ AT P n ap\^ + n ap\gap (18) 



are introduced for an arbitrary function g ap . Equations for B an ii.p) and E an i{p) are not written due to 
their complexity. The kernel K aa i{p, p') of the integral equations Jl5)— fTT) is defined by relations 

I ap [n + Sf) = £ J" d 3 p'M aa ,(p, p')8i a , p , +0{8?) , 

n ap (l + w ap ) K aa , (p, p') = -M aa > [p, p') n a , p , [l + n a , p , j . (19) 

Equations (15j— (T7) contain the values a, a defined by (8) and depending on T. The meaning of vari- 
ables v, vj, h (they are also functions of the temperature) in these equations becomes clear from the 
following expressions for the fluxes of momentum, energy, and frictional force density 

dT , 3 j 2 , 
tnl = P<>nl + Vnl,msU m U s -r] n i ms — OC n i ms U s - + 0[g A ,g A J, 

dx s dx m 

qi = hui- k^- + O (g°A 3 ,g 1 A 2 ) , Bj = -v Ul -v r ^+ if ,U + O {g° A 3 , g 1 A 2 , g 2 A 2 ) . (20) 

These expressions are obtained in the perturbation theory from the definitions (6) taking into account the 
contributions (14) to the phonon distribution function f ap [x, f ) . The term _ff| 2,1) is included in (20) because 

it gives contribution of the same order as a viscous momentum flux t^j in hydrodynamic equation for 
the velocity [see (24)] . Contributions f^, ', fj 2 ^ 1 ' to the function f ap [x,() are required in order to calculate 
b| ' (they will be analyzed in another paper). In (20) v, vj are damping rates, k is a thermal conductivity, 
Vni.ms is a viscosity tensor, h is a coefficient of the drift energy transfer, p n i ims is a coefficient of the drift 
momentum transfer, a n i itns is a coefficient of convective momentum transfer. Their values are given by 
formulas 

1 / d£ ap \ l / de ap 

h = -{Eap- A an {p)), K = --(E a p^—C an {p) 



(21) 



J \ "fn / 3 \ 

/ de ap \ _ I de ap 

Tlnl,ms — ~\ Pn~Z Dams(p) ) , Pnl,ms — \ Pn~Z Bamsip) ) , 



dpi I \ dpi 

Vn ~d~Pl £ams(p) ) ; v =^{Pn,A an (p)} U , V T = ^ Canip)}" ■ 



Here, damping rates are written in terms of a bilinear form defined by the relation 

{gap,K p } N,U = - (2 * h)3 L J d 3 pd 3 p'gapM" a "(p,p')n a , p ,[l + n a , p ,}h a , p ,, 

B 

{gap, hap} = {gap,Kp} N + {gap,Kp} U (22) 

where M^ a , (p, p'), M u , (p, p') are contributions to the kernel M aa ' (p, p') associated with the normal and 
Umklapp processes, correspondingly. Forms of this type are very useful in the kinetic theory and have 
the properties 

\gap,K p \ = {h ap ,gap\ , {gap, gap) 3=0 (23) 

that can be proved through the usual way for arbitrary functions g ap , h ap (see, for example, (U). 

Expressions (21) for v and vj follow from the integral equations (15), (16) taking into account the 
additional condition for A an {p) from (15) and expression (18) for pressure p. Formula (21) for h follows 
from equations (17) taking into account the definition of the heat capacity c (18) and the symmetry of the 



43007-5 



S.A. Sokolovsky 



bilinear form {gap, h ap } (23) . Therefore, equations (T5}-(17) can be solved without taking into account the 
expressions (21) for v,vj and h. 

The equations of phonon hydrodynamics (9) can be written in the form 



dT \ du n } dh dT lav 

C-r- - \-h— \ — ~ — ~ U n H U n 



dt { dx n J dT dx n a 



1 I dp\ dT 
" v { T dT) dx n 



dq'n 



(1,0) 



dx n 



+ 0(g°A 3 ,g 1 A 2 ,g 2 A 2 ) 



du n ( dp dT } dT dt ! .... l da dq^' 0) 

= \~ vu n-£— \-VT— -^- + R^---u n ^— + 0(g°A\g l A 2 ,g 2 A") (24) 



dt { dTdx n ) dx n dxi c dT dx\ 

with accounting relations (5), ([6), (12) , (14). Only the terms given in curly brackets were obtained, for 
example, in book [3] that discusses the standard theory [in the standard theory, h and v should be taken 
from (28)1. Other terms in (24) are obtained in the present paper (the corresponding contributions of 



the standard theory were not found in (3j). Dissipative fluxes q^'°\ t^^ are given in (20) . The obtained 



equations of the phonon hydrodynamics (24) fully take into account all the dissipative processes. 

The integral equation (15) is an equation for eigenfunction A an {p) and eigenvalue via. The solution of 
this spectral problem describes the kinetic mode announced above. According to equation (24) the velocity 
really decays if the ratio via is positive. Due to (23) , the positivity of the coefficient via follows from the 
relation 

^(Aan(p) 2 ) = {A an {p),A an (p)} (25) 

that can be derived from equation (15) and definitions (18) and (22) . This result is obtained in the frame- 
work of the considered theory which implies that the phonon drift velocity is small and has no obvious 
restrictions for the temperature. 

The considered hydrodynamic states occur when t » To where time To, introduced by the functional 
hypothesis (1). should be much smaller than the drift velocity decay time t u [t u — CT/vdue to (24)1. 



4. Comparison with the standard phonon hydrodynamics 

In the standard theory , the hydrodynamic distribution function in the zero order approximation 
in gradients ^p( x >0 coincides with the local Planck distribution with velocity n ap {^{x)). The functions 
Aanip), B an iip) according to (2), (12) . and (14) are given by expressions 

< n W = jPn, B° nl {p) = ^(l + 2n ap )pnPi. (26) 

These functions are not solutions of the integral equation (15) and the corresponding equation for B an i{p), 
i.e. A° an {p) + A an {p), B° anl {p) + B an i(p). The first inequality means that even in the linear approximation 
in drift velocity, f^p(x,f) + n ap {<;{x)) and, therefore, the local equilibrium is violated in the phonon hy- 
drodynamics. The result A° an {p) + A an (p), B° anl {p) + B an i(p) follows from the next section where it is 
shown that expressions (26) give only the leading contribution to A an {p), B an i (p) in the low temperature 
approximation. 

In the standard theory, temperature and drift velocity are defined by formulas (7) giving the following 
expressions for values a and a in (8) 

a° = ^(p 2 ), a° = -^( £ap (l+2n ap )p 2 ) (27) 

(hereinafter we denote the value A of the standard theory by A"). 

Formulas (21) can be also obtained in the standard approach but not all corresponding contributions 
to phonon hydrodynamic equations are discussed in l^-tZD - Therefore, (21) and (26) lead to expressions 
for the damping rates v, Vj and drift transfer coefficients h, [J>ni,ms 
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Equations (T6)-(T7) for the functions C an {p) and D an i{p) are simplified by the first formula in (26). 
Therefore, according to equation ( 116) and expression d2lt the phonon thermal conductivity is given in 
the standard theory by formula 

T 2 

*° = —{C°anW,C° an {p)}. (29) 



5. Solution of the integral equations of the theory at low temperatures 

At low temperatures T « Tu, the following estimates of the kernels of integral equations (15)— (17) are 
valid 

K» a ,(p,p')~p°, K u aa ,{p,p')~p\ (^e- T ° IT ) (30) 

(see, for example Hi) where Td is Debye temperature which is equal to the maximal energy of a 
phonon. 

The solutions of equations (151— (17) for functions A an (p), C an {p), andD a „;(p) and the corresponding 
equation for B an i {p), damping rates v, vj are found in the form 

A an = A M +A y n + 0(p}), v = vU +V W + 0(pi); B anl =BZ l + BM l + 0(p. 2 ); 

C an = C l ° n + C l a 1] n + 0{p. 2 ), v T = vf+vf + 0{^); D anl = Df nl + D [ l ] nl + 0{p 2 ) (31) 

(A ls] is a contribution of the order p s to the quantity A). 

It is easily to understand that the mentioned main equations of the developed theory fl"5)-(T7l are sol- 
vable by a simple iteration procedure. At each step, we obtain an integral equation of the type 

g a (p) = £ f d 3 p'K N aa , [p,p') h a , {p>) (32) 

for a function h a {p) with additional conditions that eliminate arbitrariness in h a (p) of the form c n p n + 
C£ap Igaip) is a known function]. Values v,vj, and h in this equations and kinetic coefficients of the 
system are calculated using the formulas (21) . The mentioned arbitrariness is related to conservation 
of energy and momentum in normal phonon processes. Further analysis of equations of the type (32) 
requires additional information on the spectrum e ap of the phonons and amplitude of their interaction 
<$(.aip\;a2P2,o:3P3). It should be also stressed that necessary solutions of equation (32) are vectors and 
tensors of second rank depending on momentum. Even for crystals of the considered cubic symmetry 
classes O, Oh, Tfc, they have a rather complicated structure. 

Note also that at low temperatures, the long-wavelength acoustic phonons give the main contribu- 
tion to thermodynamic and kinetic properties of a crystal (see, for example, [6jj). This allows one to find 
phonon spectrum e ap and phonon interaction amplitude <£>{a\p\\a2P2, 0:3^3) based on the elasticity the- 
ory. This leads to some simplification of the above obtained equations that will be analyzed elsewhere 
(the standard theory in this approach has been constructed in 

Let us present some results of our calculations in the main low temperature approximation: 

AZW = -^Pn, v™ = -2—{ Pn , Pn } v , H m = 

(p 2 ) (p 2 ) (p 2 ) dT 

V™ = ±{Pn,C2( P )} U , *™ = ^{Cf n ip),C^p)f, 



(P 2 ) 



'nl,ms 3 fj 



{Df nl ( P ),Df ms {p)} N . (33) 



Among other things, these formulas give positively defined expressions for the heat conductivity and 
viscosity of the system. Some results of the main approximation coincide with the corresponding results 
of the standard theory 
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v [i] =v o h M =h O t ^l ms = ^ ims> (ata.aV^A (34) 

Functions C% n , D° an , satisfy the integral equations (16), (T71 with quantities A° an , a", h°,v°, v° taken from 
(26l-(28) instead of A an ,a, h,v,vj. These equations can be solved at low temperatures in the perturba- 
tion theory in p that gives the following results 

cZw-cTw, D>=flffl(P)i 

4 ] =xf\ k [01 =k° |01 , V ] =T}°f , (at a = a°, a=a°) (35) 
1 1 ' 'nl,ms 'nl,ms' v 

(notation like i4 0|s| gives a contribution of the order p s to a quantity of the standard theory A ). So, results 
of the developed theory and the standard theory coincide with one another in the main low temperature 
approximation (at o - a", a = a , i.e., at the standard definition of the drift velocity and temperature). 

The first two formulas (34) show that at low temperatures, the difference between the phonon dis- 
tribution function in the zeroth approximation in gradients f™ [x, <f) and the local Planck distribution 
n a p{£,{x)) is estimated by f£p(je,f) = n ap {^{x)) + 0{A, 1 p l , A 2 p 1 , A 3 ). Therefore, at low temperatures T« To 
and small drift velocity, the violation of the local equilibrium in the phonon hydrodynamics is exponentially 
small. This result confirms the applicability of the standard theory at low temperatures. 



6. The steady states of insulators 

Investigation of steady states of the system is of great interest because their properties are easier to 
be analyzed experimentally. The hydrodynamic equation for the drift velocity (24) should be solved in 
the steady state with respect to the drift velocity u n in the form of a series in temperature gradients that 
gives 

1 ( dp) dT n( 3 



vh + ^J^ + ^- (36) 



The accuracy of this result is limited by the accuracy estimation in (24). It should be noted that an equation 
of the type (36) appears in the standard theory of steady processes in dielectrics too [see, for example, 0] 
where (36) is obtained for T « Tq without a term with vj that is exponentially small according to (31)1. 
Taking into account the result (36) . the expression (20) for phonon energy flux yields 



«" = -*S +0(g3) ' *" + £( vr+ !f) (37) 

where we have introduced thermal conductivity in a steady state ic. The formula shows that in the ab- 
sence of the Umklapp processes, conductivity k = oo [according to (21) , in this case v = 0, vj — 0]. This 
well-known result is not surprising since an isolated system does not have nonequilibrium steady states. 
Only the presence of Umklapp processes that lead to a nonconservation of the phonon momentum of a di- 
electric crystal makes steady states possible. 

There is an interesting issue concerning the conditions of stationarity of the phonon temperature of 
such a system. Substitution of the expression for the drift velocity in a steady state (36) into equation (24) 
for temperature, yields a condition 

d LdT\_ Q 
dx„ \ dx n I 

(in the second order of the perturbation theory in gradients). The meaning of this result as the condition 
of heat source absence in a steady state is clear. At the same time, the issue regarding the temperature 
distribution in isolated insulators is not discussed in the literature, although the nature of steady states 
in such systems is quite unusual. 

Steady state thermal conductivity of the phonon system of an insulator k can be calculated at low 
temperatures based on the results of the previous section that in the main approximation gives 

* = *^ + 0(A K^ = -^-m 2 . (39) 

\Pn,Pn\ y0l > 
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This result coincides with the result by Akhiezer 0] and does not depend on the definition of the drift veloc- 
ity (8). The thermal conductivity k is exponentially large which eliminates the problem of calculating 
the corrections to it at T « Tv. 

7. Conclusions 

In this paper, the Chapman-Enskog method is generalized to take into account relaxation processes 
(kinetic modes) in the hydrodynamic theory, i.e., the processes that can be present in a spatially homoge- 
neous state of a system. On this basis: 

• Nonlinear hydrodynamics of the phonon system has been built in the approximation of small 
phonon drift velocity for insulators with cubic symmetry of the lattice . 

• It is proved that a small phonon drift velocity decays. 

• In the perturbation theory in the drift velocity it is found, that in the phonon hydrodynamics the 
distribution function of phonons in the zeroth approximation in the gradients is different from the 
Planck distribution with the velocity. At the same time, at zero drift velocity, these distributions 
coincide. So, it is established that the local equilibrium in the phonon hydrodynamics is violated. 

• It is shown, that at low temperatures the integral equations of the theory are solvable iteratively. 

• It is shown, that at low temperatures and at small drift velocity, the local equilibrium in the hydro- 
dynamic phonon system takes place with exponential accuracy. 

• The obtained results can be applied to the analysis of hydrodynamic processes in the system of 
phonons of the insulator at intermediate and high temperatures. 

• The relation between usual hydrodynamic thermal conductivity and thermal conductivity of phonons 
in the steady state has been established. It is shown that the Akhiezer expression gives the main 
low temperature contribution to the thermal conductivity of an insulator in its steady states. 

• The thermal conductivity of an insulator is impossible to be calculated without taking into account 
the Umklapp processes because nonequilibrium steady states in a closed system do not exist with- 
out them. 
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riflpoflMHaiviiMHi cTaHi/i <|>OHOHiB AiesiexTpi/iKiB 

CO. COKOTlOBCbKI/M 

npuflHinpoBCbKa flepxaBHa aKafleiviifl 6yfliBHWLrrBa Ta apxiTeioypn 

MeTOfl MenMeHa-EHCKOra y3ara/ibHeHO pj\n BpaxyBaHHfl Bn/iwBy KiHeTUHHUx MOfl cucTeMW Ha riflpoflHHaMHHy 

eBO/1K)L4iK). y Ha6/lWXeHHi Ma/lOi flpefi(|>OBOi LUBWflKOCTi BUBMeHO TiflpOflHHaMHHi CTaHU (flOHOHHOi niflCWCTe- 

mu flie/ieicrpuKa. fl,nn cnpoLueHHfl, floaiiflxeHHfl npoBefleHO fl/ia KpucTa/iiB Ky6i4Hnx K/iaciB cwMeTpii. Bcra- 
HOB/ieHO, mo y cfiOHOHHiw riflpoflWHaMiL4i /lOKa/ibHa piBHOBara nopywyeTbCfl HaBiTb y /liHMHOMy Ha6/iwxeHHi 
3a LUBUflKk™. I4e e Hac^iflKOM BiflcyTHOCTi 3ai<OHy 36epexeHHfl iMny/ibcy, mo Befle flo pe/iaKcaijiT flpeRc|>OBoT 
LUBHflKOCTi. OTpnMaHO piBH^HHfl <j>OHOHHoi riflpoflHHaMiKW 3 ypaxyBaHHHM flncwnaTnBHUx npoijeciB. Pe3y/ib- 
Ta™ nopiBH^HO 3i CTaHflapTHOK) Teopiero, flKa 6a3yeTbcn Ha HanBHOcri /lOKa/ibHOi piBHOBa™. OTpuMaHO iHTe- 

rpa/lbHi piBHflHHfl fl/lfl p03paxyHKy 06'EKTiB Teopii (BK/lKDHaK)HH B'fl3KiCTb Ta Ten/lOnpOBiflHiCTb). noKa3aHO, luo y 

rpaHHL4i HH3bKnx TeMnepaTyp uj piBHAHHfl po3B'H3yKnbCH iTepai_(iflMH. P03i7iHHyT0 CTai^iOHapHi CTaHU cucTeMU 
i OTpnMaHO Bwpa3 fl/ia IT CTai4iOHapHoT TenflonpoBiflHOdi. noKa3aHO, luo BOHa y HH3bKOTeMnepaTypHifi rpaHWL(i 
cniBnaflae 3 BiflOMHM pe3y/ibTaTOM Axie3epa. BcTaHOB/ieHO, luo po3nofli/i TeMnepaTyp y CTai4ioHapHOMy CTam 
flie^eKTpwKa 3aflOBO/ibHae yMOBy BiflcyTHOCTi flxepe/ia Ten/ia. 

IOi K) h o b i c/iOBa: tyOHOHW AieneKTpvma, npou,ecvi nepem^anm, penaKcaujviHt cryneHi Bi/ibHOcri, noKanbna 
piBHOBara, Merofl MenMeHa-EHCKOra, cpOHOHHa riflpoflHHaMiKa, Ma/ia ppevupOBa LUB^KiCTb, HW3bK\ 
rewmeparypw, craLfiOHapHi cranm 
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